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Abstract 

This paper investigates the achievable sum rate of multiple-input multiple-output (MIMO) wireless systems 
employing linear minimum mean-squared error (MMSE) receivers. We present a new analytic framework which 
unveils an interesting connection between the achievable sum rate with MMSE receivers and the ergodic mutual 
information achieved with optimal receivers. This simple but powerful result enables the vast prior literature on 
ergodic MIMO mutual information to be directly applied to the analysis of MMSE receivers. The framework is 
particularized to various Rayleigh and Rician channel scenarios to yield new exact closed-form expressions for the 
achievable sum rate, as well as simplified expressions in the asymptotic regimes of high and low signal to noise 
ratios. These expressions lead to the discovery of key insights into the performance of MIMO MMSE receivers 
under practical channel conditions. 
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I. Introduction 

Multiple-input multiple-output (MIMO) antenna wireless communication systems have received enor- 
mous attention in recent years due to their ability for providing linear capacity growth without the need 
for increased power and bandwidth [1,2]. Since the important discoveries in [1,2], a major focus has 
been directed at investigating the MIMO channel capacity under a wide range of propagation scenarios. 
For example, the impact of physical phenomena such as spatial correlation, line-of-sight, antenna mutual 
coupling, frequency-selectivity, and co-channel interference, have now been well-studied, especially for 
single-user MIMO systems [3-18]. Many key results have also been derived in the multi-user context, 
and this is still an important topic of on-going research (see, eg. [19-23], and references therein). 

Despite the abundance of literature on MIMO channel capacity, the vast majority of existing work in 
this area has focused primarily on systems employing optimal nonlinear receiver structures. It is very well- 
known, however, that such receivers can have prohibitive complexity requirements for practical systems 
and that low complexity receivers, such as linear minimum mean-squared error (MMSE) receivers, are 
more appealing for many applications. 

However, despite their practical importance, there are currently few closed-form analytical results on the 
achievable rates of MIMO MMSE receivers. In fact, most prior related work has focused on studying the 
asymptotic achievable rates of linear MMSE receivers, and has been derived in the context of code-division 
multiple access (CDMA) systems employing random spreading. It is well-known that such systems, under 
certain conditions, are isomorphic to single-user MIMO MMSE systems. In this context, the primary 
approach has been to study the asymptotic spectral efficiency and signal to interference plus noise ratio 
(SINR) as the system dimensions grow large with fixed ratio, using advanced tools from large-dimensional 
random matrix theory [24-29]. Those results, which are extremely accurate for the high system dimensions 
encountered in CDMA applications (determined by the spreading factor and the number of users), may 
become less accurate for the much smaller system dimensions indicative of current MIMO systems 
(determined by the numbers of transmit and receive antennas). Moreover, in many cases, the existing 
CDMA results restrict the equivalent channel gains to be independent random variables with zero mean; 
thereby precluding many MIMO statistical channel models of practical interest (eg. correlated Rayleigh 
and Rician fading). An exception is the very recent work [30], which investigated the asymptotic (large 
antenna) mutual information distribution of MIMO MMSE receivers in the presence of spatial correlation. 



For finite-dimensional systems with arbitrary numbers of transmit and receive antennas, there are 
remarkably few analytic results dealing with the achievable rates of MIMO MMSE receivers. The only 
directly related results, of which we are aware, were derived very recently in [23, 31, 32]. Specifically, [31] 
investigated the achievable rates of MIMO MMSE receivers by characterizing the asymptotic diversity- 
multiplexing trade-off. In [32], expressions for the exact achievable sum rates were derived for MIMO 
MMSE receivers for uncorrelated Rayleigh fading channels, based on utilizing the distribution of the 
corresponding SINR at the MMSE receiver output. A similar method was employed in [23], which 
presented analytic expressions for the high signal to noise ratio (SNR) regimqj, again restricting attention 
to uncorrelated Rayleigh channels. The derivation approaches in [32] and [23] appear intractable for more 
general channel models. In [33], properties of mutual information were used to conclude that the MMSE 
receiver with perfect decision feedback (ie. nonlinear) is optimal in certain scenarios. However they did 
not consider the realistic case of non-perfect feedback, nor the basic case of no-feedback (ie. the linear 
MMSE receiver); and they did not consider fading. 

In this paper, we introduce a new general analytic framework for investigating the achievable sum 
rates of MIMO systems with linear MMSE receivers. Our main results are based on some very simple 
but extremely useful algebraic manipulations which essentially relate the MIMO MMSE achievable sum 
rate to the ergodic MIMO mutual information with optimal receivers. This relationship permits us to 
circumvent the extreme difficulties entailed with explicitly characterizing the SINR distribution at the 
output of the linear MMSE filter, and instead to directly draw upon the vast body of existing results on 
ergodic mutual information from the MIMO literature. In particular, using this general framework, we 
can obtain analytic expressions for the achievable sum rate of MIMO MMSE receivers in a broad class 
of channel scenarios, without the need to invoke a large number of antennas. 

We demonstrate our approach by first considering the canonical uncorrelated Rayleigh fading channel, 
which, as mentioned above, has already been tackled using different methods in [32] and [23]. We show 
that by employing our new framework, it is easy to obtain equivalent expressions to those presented in 
[32] and [23], in addition to establishing new results. We then consider spatially-correlated Rayleigh and 
uncorrelated Rician (line-of-sight) fading channels, for which there are no comparable prior results. For 

'Note that [23] considered the different context of multi-user MIMO broadcast channels with linear zero-forcing precoding; however, for 
the high SNR regime, there is a strong analogy to the single-user MIMO MMSE model considered in this paper. 



these channels, by employing our general analytic framework, we derive new exact expressions for the 
achievable sum rates of MIMO MMSE receivers, as well as simplified characterizations in the high and 
low SNR regimes. 

In many cases, our new analytical expressions show a clear decoupling of the effects of transmit 
correlation, receive correlation, and line-of- sight, which leads to key insights into the performance of 
MIMO MMSE receivers under practical channel conditions. For example, they reveal the interesting 
result that at both high and low SNR, the achievable sum rate of MIMO MMSE receivers is reduced 
by either spatial correlation or line-of-sight. At high SNR, this rate loss is due to an increased power 
offset, whereas at low SNR, through a reduced wideband slope. We also present an analytical comparison 
between the achievable rates of MIMO MMSE receivers, and those of optimal receivers. Interestingly, we 
show that at both high and low SNRs, although both MMSE and optimal receivers incur a rate loss due 
to either spatial correlation or line-of-sight, the loss is more significant for MMSE. 

The paper is organized as follows. In Section HI] we introduce the basic MIMO MMSE signal model of 
interest, and its corresponding achievable rate. We then present our general analytic framework in Section 
HIH before particularizing these general results to various Rayleigh and Rician fading models in Section 
HVl Finally, Section |V] gives some concluding remarks. 

II. Signal Model and Achievable Sum Rate 

Consider a single-user MIMO system with A^^ transmit and A^^ receive antennas, with discrete-time 
input-output relation 

r = Ha+n (1) 

where r is the A^^ x 1 received signal vector, n is the A^^ x 1 vector of additive white Gaussian noise with 
covariance E [nnT] = A'oIat^, and a is the A^t x 1 vector of transmit symbols, satisfying the total power 
constraint E [a^a] = P. The A^ x Nt matrix H represents the flat-fadingn MIMO channel, assumed to 
be known perfectly at the receiver but unknown to the transmitter, and normalized to satisfy 

£;H[tr(HHt)] =NrNt. (2) 

^If the fading is frequency-selective, then our results can also be easily applied upon decomposing the channel into a set of parallel 
non-interacting frequency-flat subchannels. 



Throughout the paper, we consider the class of MIMO spatial multiplexing systems with independent 
equal-power Gaussian signalling, in which case the input signals have covariance E [aa^] = -^iNf 
For optimal receivers, the ergodic mutual informatiory is given by 

/°P*(snr, Nr, Nt) = En [P^\sm, N,, N^, H)] (3) 

where snr = P/Nq, and 

+ — HHM . (4) 

In this paper, we focus on characterizing the achievable sum rate of linear MMSE receivers. Such 
receivers operate by applying a linear filter to the received signals to form the estimate 



a ^^mmsel" ^^mmse (Ha -|- II) 



with Wjnmse choscu to minimize the mean-square error cost function 

Wmmse = aig Hiln E [||a- Grf] . 
The solution to this optimization problem is well-known (see, eg. [34]), and is given by 



W 



P 
P 



HHt + ^I^^ 
snr 




HtH + ^I^, 



snr 



-1 



Ht 



(5) 

(6) 

(7) 
(8) 



where the second line is due to the matrix inversion lemma. It can be easily shown (see eg. [35]) that the 
instantaneous received SINR for the ith filter output (ie. corresponding to the ith element of a) is given 
by 

1 



7i 



1 



(9) 



with [-jj j denoting the ith diagonal element. Assuming that each filter output is decoded independently, 
the achievable sum rate is expressed as 



Nt 



J— (snr, iV,, N,) = ^ E,^ [log, (1 + 7.)] . 



(10) 



i=l 



In general, the exact distribution of 7^ does not appear to be available in closed-form, other than for 
the specific cases of independent and identically distributed (i.i.d.) or semi-correlateqj Rayleigh fading 

^Note that for the case of uncorrelated Rayleigh fading channels, this corresponds to the ergodic capacity. 
''The term semi-correlated refers to channels with correlation at either the transmitter or receiver, but not both. 
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[36]. This precludes direct evaluation of (flOl) for many channels of practical interest. 

III. General Analytic Framework for the Achievable Sum Rate of MIMO MMSE 

Receivers 

In this section, we present our new general analytical framework for investigating the achievable sum 
rates of MIMO systems with linear MMSE receivers. In particular, we show that by using some very 
simple manipulations, the achievable sum rate (fTOl) can be expressed in a form which can be easily 
evaluated for many fading models of interest, without requiring explicit statistical characterization of 7^. 

We find it convenient to introduce the following notation: n = mm{Nr, Nt), m = max{Nr, Nt), n' = 
min(A^^, Nt - 1), and m' = max(A^,., Nt-1). 



A. Exact Characterization 

The following key result presents a simple and powerful connection between the MIMO MMSE 
achievable sum rate and the ergodic mutual information obtained with optimal receivers. 



(11) 



Theorem 1: The achievable sum rate of MIMO MMSE receivers can be expressed as 

^* r /TV - 1 

I~isnT,Nr,Nt) = NtEn [/°P*(snr, iV„ iV*, H)] - J^^h, 1°"' f ^^snr, iV„ iV, - 1,H 

where Hj corresponds to H with the ith column removed. 

Proof: See Appendix HI ■ 
Corollary 1: When H contains i.i.d. entries, ([TT)) reduces to 

J— (snr, Nr, Nt) = Nt (^/"^^(snr, iV„ Nt) - P^' (^^^^nr, iV„ iV, - 1^ ^ . (12) 
Importantly, with the MMSE achievable sum rate expressed in this form, the required expectations 
are the same as those required for the evaluation of the ergodic MIMO mutual information with optimal 
receivers which, as already discussed, have well-known solutions for many channels of interest. In Section 
HVlwe will draw upon these previous results to yield new closed- form expressions for the MMSE achievable 
sum rate. 



B. High SNR Characterization 

In the high SNR regime, the ergodic MIMO mutual information and the achievable sum rate of MIMO 
MMSE receivers can be expressed according to the affine expansioijfl [25] 

J(snr, Nr, Nt) = Soo (logs snr - C^o) + o(l) (13) 

where Soo is the high SNR slope, in bit/s/Hz/(3 dB) units, given by 

I(snT,Nr,Nt) 

Soo= lim \ ' ^' (14) 

snr^oo log2 Sni 

and Coo is the high SNR power offset, in 3 dB units, given by 

/ Iism,Nr,Nt) \ 
Coo= lim (logasnr 1 . (15) 



snr— >oo 



For MIMO systems with optimal receivers, these parameters are obtained from ([3]) as 

ST = mHNr,Nt) (16) 

and 



= log2 Nt - -En [J{Nr, Nt, H)] (17) 

1% 



respectively, where 



logo det (HHt) , Nr<Nt 
J{N^,Nt,U) = { ^ ' ' . (18) 

log2 det (HtH) , Nr>Nt 
For MIMO systems with MMSE receivers, we have the following key result: 

Theorem 2: At high SNR, the achievable sum rate of MIMO MMSE receivers can be expressed in the 
general form (fT3l) with parameters 

f Nt, for Nr > Nt 
5— =<^ (19) 
[ 0, for Nr < Nt 

and C^^^ is given by 

~ _ i Nt - En [J{Nr, Nt, H)] + ^ E£i ^h, [J{Nr, Nt-1, H,)] , N, > Nt 

— \ ■ (2U) 

[ OO, Nr < Nt 

Proof: The result is easily obtained by substituting (fTTI) into (fT4l) and (fT5l) and evaluating the 
necessary limits. ■ 

^The notation f{x) — o(g(x)) implies tliat lim^i^cx) ~ 0. 



Corollary 2: Consider the case A^^ > Nf. If H has i.i.d. entries, then (|20l ) reduces to 

C~ = log2 Nt - Eu [J{Nr, Nt, H)] + Eu, [J{Nr, Nt - 1, Hi)] . (21) 
The fact that 5™™'''^ = and = oo for the case A^^ < A^^ occurs since, when decoding the 

data stream for a given transmit antenna, the MMSE receiver does not have the capabilities (or enough 
degrees of freedom) to perfectly cancel the interference caused by the other transmit antennas. Thus, even 
when the impact of receiver noise becomes negligible (i.e. as p ^ oo), the channel remains interference- 
limited and, as expected, the MMSE achievable sum rate converges to a finite asymptote. For the more 
interesting case, with > Nt, we see that the high SNR power offset is non-zero. Very importantly, 
this case involves the same types of expectations as those required for the high SNR analysis of ergodic 
MIMO mutual information, for which closed-form solutions exist for a wide range of fading channel 
scenarios [5, 15,37,38]. In the sequel, we will draw upon these previous results in order to derive new 
simple closed-form expressions for under a range of conditions. 

In addition to the absolute high SNR power offset, it is also of interest to to examine the excess high 
SNR power offset with respect to the ergodic MIMO mutual information Considering the case 

Nr > Nt, this is given by 

Aex = C~ - CT ■ (22) 

This measure is meaningful, since both the MMSE and optimal receivers yield the same high SNR 
slope (under the assumption that Nr > Nt), and, as such, the corresponding curves will be parallel. 

C. Low SNR Characterization 

When considering the low SNR regime, it is convenient to introduce the concept of the dispersion of a 
random matrix. This measure, originally introduced in [7], will play a key role in subsequent derivations. 
Definition 1: Let © denote a N x N random matrix. Then the dispersion of is defined as 

CO) = iV^M^ . (23) 
For low SNR, it is often appropriate to consider the achievable rate in terms of the normalized transmit 
energy per information bit, rather than per-symbol SNR. This can be obtained from /(snr) via 



J ( 1^ ) = /(snr) (24) 
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with snr the solution to 

E}, snr 



A^o /(snr) 

Note that ^ is related to the normalized received energy per information bit, via 



(25) 



In general, closed-form analytic formulae for (|24l) are not forthcoming, however, for low levels, this 
representation is well approximated by [27] 

xf#l-^olog,f^l (27) 



No 



mm - 



where the approximation sharpens as ^ [ ^^j^- Here, ^^^^ and S^, are the two key parameters which 
dictate the behavior in the low SNR regime corresponding, respectively, to the minimum normalized energy 
per information bit required to convey any positive rate reliably, and the wideband slope. Importantly, 
they can be calculated directly from /(snr) via [27] 

Eb snr 
— — = lim 



^Omin s°=^-^o /(snr) 

^ (28) 



m 

and 

So= lim — p ^ °\ ^101ogio2 

^ I ^ 10 logio - 10 logio 

= in2 (29) 

1(0) 

respectively, where /(•) and /(■) denote the first and second-order derivative respectively, taken with 
respect to snr. Note that I (^jf^^ implicity captures the second-order behavior of /(snr) as snr 0. For 
MIMO systems with optimal receivers, these parameters are given by [7] 

= (30) 

' mm ' r 

and 

2K 



So"' = TTTT^ (31) 



C(HHt) 
respectively. 

For MIMO systems with MMSE receivers, we have the following key result: 
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Theorem 3: At low SNR, the achievable sum rate of MIMO MMSE receivers can be expressed in the 
general form (|T7l) with parameters 

In 2 

. = AT (32) 



and 



S-Q™ = '^^^ . (33) 

A^,C(HHt) - (^) Ef:iC(H.Hl) 



Proof: See Appendix [III 
Corollary 3: If H has i.i.d. entries, then (l33l) reduces to 



2A^ 

° " iV,C(HHt) - li^sFCCHiHl) ■ ^ 

Interestingly, comparison of (l30l) and (|32|) reveals that MMSE receivers are optimal in terms of the 
minimum required For both receivers, this parameter is independent of the number of transmit antennas, 
whilst varying inversely with the number of receive antennas; a fact directly attributed to the increased 
channel energy captured by the additional receive antennas, whilst the total transmit energy is constrained. 
We also see that the wideband slope of MIMO MMSE receivers depends on the random matrix channel 
via its dispersion. In the following section we will evaluate this parameter in closed-form for various 
Rayleigh and Rician fading channels of interest. From these results, we will see that although MMSE 
receivers are optimal in terms of the minimum required such receivers are indeed suboptimal in the 
low SNR regime as typically reflected in a reduced wideband slope 5*0. 

IV. Achievable Sum Rate of MIMO MMSE in Fading Channels 

In this section, we demonstrate the key utility of the general results propounded in the previous section, 
by presenting explicit solutions for the MIMO MMSE achievable sum rate for various fading models of 
practical interest. 

A. Uncorrelated Rayleigh Fading 

We start by particularizing the results for the canonical case: the i.i.d. Rayleigh fading channel, 

H ~ CMn^,n, (0, In^ ®In,) , (35) 
representative of rich scattering non-line-of- sight environments with sufficiently spaced antenna elements. 
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SNR (dB) 



Fig. 1. Achievable sum rate of MIMO MMSE receivers in i.i.d. Rayleigh fading; comparison of analysis and simulations. Results are shown 
for different Nr = Nt = n. 

1 ) Exact Analysis: 

Proposition 1: For i.i.d. Rayleigh fading, the MIMO MMSE achievable sum rate is given by 

I {snr,N.,N,)-N,e y r„(m)r„(n) TAm')TAn') j ' ' ^^^^ 

where r„(-) is the normalized complex multivariate gamma function, 

n 

TnM = J]r(m-^ + l) (37) 

1=1 

and ^n.mik) is an n x n matrix with (s,t)th element 

f ^.t'EI-r'Ej^) for t = fc 

{^n,mm,,= { ^'""'^ (38) 

[ Ts,t^. for k 

where Ts^t = n + m — s — t, and E/i(-) is the exponential integral. 

Proof: This result is easily obtained by_plugging into (fT2l) the ergodic mutual information expression 

[39]. ■ 



for i.i.d. Rayleigh MIMO channels given ire 



6 



An alternative closed-form expression for ergodic mutual information can be found in [5]. 
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Our result in Proposition \T\ gives an exact closed-form expression for the MMSE achievable sum rate, 
which applies for all SNRs and arbitrary antenna configurations. This result is confirmed in Fig. [H where 
it is compared with the exact MMSE achievable sum rate, obtained via Monte-Carlo simulations, for 
different antenna configurations. There is precise agreement between the simulated and analytic curves, 
as expected. We note that Proposition\l}pr:esents a new expression for the achievable sum rate of MIMO 
MMSE receivers, however, an alternative expression has also been obtained via different means in [32]. 
That result was obtained by directly integrating (flOl) over the distribution of the SINR in Q; an approach 
that cannot be followed for more general channel models. 

For small system dimensions (eg. n = 2), (l36l) reduces to particularly simple forms. For example, for 
the case Nt = 2, Nr > 2, it reduces to 



/~(snr, Nr, 2) = 2e'/'^' ^ E, (2/snr) + iV, (E^,+i (2/snr) - E^, (2/snr)) log^ e , (39) 



2) High SNR Analysis: Recall that in the high SNR regime, the key channel-dependant parameter is 
the high SNR power offset for the case Nr > Nt. 

Proposition 2: For i.i.d. Rayleigh faded channels, the high SNR power offset (for A^,. > A^^^) for MIMO 
MMSE receivers is given by 




k=i 



whilst for Nr = 2, Nt > 2, we get 




(40) 




(41) 



where 7 ^ 0.5772 is the Euler-Mascheroni constant. 



For Nr = Nt 



t — 



n this reduces to 



mmse 



log2(rae'^) 



(42) 



'00 



Proof: The result is easily obtained from (f2T)) . upon noting that [37] 



Nt-l 




(43) 



£=0 
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for Nr > Nt, where is the digamma function defined as 

/ for J>1 

iP{j) = ( . (44) 

[ -7 for j = 1 

U 

Recalling that the MMSE receiver behaves equivalently to the linear zero forcing (ZF) receiver at high 
SNR, we note that Proposition |2] could also be easily derived by starting with the high SNR MIMO ZF 
sum rate expression presented in [40, Eq. 8.54] for the case of i.i.d. Rayleigh fading channels. 

Together with (fT9l) . Proposition |2] indicates that if the number of transmit antennas is kept fixed and 
the number of receive antennas are increased, then, whilst having no effect on the high SNR slope, the 
high SNR achievable rate is improved through a reduction in the power offset. Intuitively, this is due 
to the additional received power captured by the extra antennas, and also to the enhanced interference 
cancelation capabilities afforded by the extra degrees of freedom in the receive array. In fact, as A^,. ^ Nt, 
then C^^^ — > — oo dB, confirming the intuition that the MMSE receiver completely mitigates the effect 
of fading (in the high SNR regime) as the number of degrees of freedom at the receiver greatly exceed 
the number of impeding interferers. 

It is also worth noting that, based on (|4T)) and (|42)) . one may conclude that increasing the number of 
transmit and receive antennas, whilst keeping their difference fixed, may have a deleterious effect on the 
achievable rate due to an increased high SNR power offset; especially when A^,. = Nt = n. However, care 
must be taken when interpreting this result. In particular, since the high SNR slope (fT9l ) also increases 
linearly with A^^, it turns out that the overall MMSE achievable sum rate actually increases with n. This 
result is seen in Fig. [H where the high SNR MMSE achievable sum rate based on ((42l) and (fT3l) is 
presented for n = 2 and n = A. We see that the slope is greatest for the case n = 4, as expected; however 
the power offset, which determines the point at which the high SNR linear approximation intersects with 
the horizontal SNR axis, is smallest for the case n = 2. 

As an aside, it is also important to note that although the general approximation (fT3]) is formally valid 
in the regime of very high SNRs, Fig. [T] demonstrates good accuracy even for moderate SNR values (eg. 
within 20 dB). 

As the next result shows, the high SNR power offset (|4T]) admits further simplifications in the "large- 
antenna" regime. 
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Corollary 4: For i.i.d. Rayleigh faded MIMO channels, as the number of antennas grows with ratio 
P = ^ (with (3 < 1), the high SNR power offset dH]) converges to 



'mmse 



log2 ( ) • (45) 



1-P 

Proof: The result is easily established upon noting that 

n— 1 

XI 7 ~ ^ = ^('^) ~ l^(^) ("^6) 

£=1 ^ 

for large n. ■ 
Interestingly, we see that the high SNR power offset is unbounded for /? = 1 (i.e. A^^ = Nt); however, 
it converges for all /3 < 1, decreasing monotonically in p. We note that this expression agrees with a 
previous large-system result derived for MMSE receivers in the context of CDMA systems with random 
spreading [25]. 

Corollary 5: For i.i.d. Rayleigh faded MIMO channels, the excess high SNR power offset is given by 

Aex = log,eU^ E --l). (47) 

\ e=Nr-Nt + l / 

For Nr = Nt = n, this reduces to 

Aex = log2e^^. (48) 
Proof: The result is obtained by substituting (|4TI) and [38, Eq.l5] 

\ e=i * e=Nr-Nt+i 



(49) 



into (|22|) . and performing some basic algebraic manipulations. ■ 
Note that an alternative expression for (l47l) can also be obtained from [23, Theorem 2] and [23, Eq. 

(15)], which considered the asymptotic excess rate offset of linear precoding in uncorrelated Rayleigh 

fading MIMO broadcast channels. 

The excess high SNR power offset also admits a simplified characterization in the large-antenna regime. 



Corollary 6: For i.i.d. Rayleigh faded MIMO channels, as the number of antennas grows with ratio 
P = ji- (with P < 1), the excess high SNR power offset (|47l) converges to 

Acx ^ log2 - log2 e . (50) 

Proof: The result is derived trivially from (1471) upon employing (|46l) . ■ 
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Again, we note that this expression agrees with a previous large-system result derived in [25], which 
considered the context of randomly- spread CDMA systems. 

3) Low SNR Analysis: Recall that in the low SNR regime, the key channel-dependant parameter is the 
wideband slope. 

Proposition 3: For i.i.d. Rayleigh faded channels, the wideband slope for MIMO MMSE receivers is 
given by 

^0 - 2Nt + N,.-l • ^^^^ 
Proof: For i.i.d. Rayleigh fading, using [7, Lemma 6], we find that 

C (HH.) ^ C(H.Hl) ^ ^^^±^ . (52) 



Substituting ([521) into dMl) leads to the result. ■ 
This agrees with a recent result obtained via different methods in [32]. It is interesting to compare (|5T|) 
with the corresponding wideband slope for optimal MIMO reception, given in [7] as 

- ■ ^^^^ 

In Fig. |2l the low SNR achievable rate approximations for MMSE and optimal receivers are presented. 



based on (|5T]) and (1531) respectively. The curves are shown as a function of received for a system 
with Nr = Nt = 3. In both cases, the corresponding exact low SNR curves are also presented for further 



comparison, obtained by numerically solving (1241) and (125]) . The figure shows that the linear approximations 
are accurate over a quite moderate range of values, especially for the MMSE receiver. 

Clearly, from (ISTl) . 5*™™^° is increasing in both Nr and Nt, with the rate of increase being more 
significant for A^^- This is in contrast to Sq^^, in which case both A^^ and Nt play symmetric roles. We 
also see that 

^mmse ^ Nt + Nr 



2Nt + Nr-1 

This ratio is increasing in A^,, and decreasing in Nt, satisfying 



(54) 



1 Nt + Nr 

- < — < I (55) 

2 - 2Nt + Nr-l - 

where the lower bound is approached as A^t oo for fixed A^^^ confirming that, relative to optimal 
receivers, MMSE receivers take a significant hit in the low SNR regime when the number of transmit 
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Fig. 2. Comparison of the spectral efficiency of a MIMO system with optimal and MMSE receivers in i.i.d. Rayleigh fading channels. 
Results are shown as a function of received for Nr = Nt = 3. 



antennas exceed the number of receive antennas. This is due primarily to the limited interference suppres- 
sion capabilities of the receive array in this "overloaded" scenario. On the other hand, the upper bound is 
achieved strictly for A^^^ = 1. It is also approached as A^^ ^ oo for fixed Nt, revealing the intuitive notion 
that linear MIMO MMSE receivers perform near-optimally if the number of receive antennas are much 
larger than the number of transmit antennas, due, once again, to the additional captured received power 
and the enhanced interference suppression capabilities of the receive array. 
Finally, it is interesting to consider the large-antenna regime. 

Corollary 7: For the i.i.d. Rayleigh faded channel, as the number of antennas grows with ratio /? = 



the ratio between the MMSE wideband slope (1511) and the optimal wideband slope (1531) converges to 



Cmmsc 
•^0 



1 + 2/3 ' 



which, interestingly, for /3 = 1 (i.e. A^^ = Nf) gives 



Cmmse 

QOpt 



2 
3 



(56) 



(57) 
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B. Correlated Rayleigh Fading 

We now particularize the general results of Section [ni] to spatially-correlated Rayleigh fading channels, 
representative of non-line-of- sight environments with a lack of scattering around the transmitter and/or 
receiver, or with closely spaced antennas (with respect to the wavelength of the signal). We consider the 
popular "separable" correlation model, described by 

^r^CMN.^N, (0,R® S) , (58) 

where R and S are Hermitian positive-definite matrices which represent, respectively, the receive and 
transmit spatial correlation. This model, commonly adopted due to its analytic tractability, has also been 
confirmed through various measurement campaignj^ (see, eg. [42,43]). 

It is important to note that, to our knowledge, the results in this section present the first analytical 
investigation of the achievable sum rate of MIMO MMSE receivers in the presence of spatial correlation. 

1) Exact Analysis: For our exact analysis, we focus on semi-correlated scenarios, allowing for spatial 
correlation at either the transmitter or receiver (but not both). We note, however, that the same approach 
can also be applied to derive closed-form solutions for the more general case in (l58l) . i.e. allowing 
for correlation at both the transmitter and receiver, by employing the MIMO ergodic mutual information 
results for such channels established in [11,44]. The final expressions, however, involve more cumbersome 
notation compared with the semi-correlated results, and as such we choose to omit them here. (Note that 
the more general model (1581) will be explicitly considered in the following subsections, when focusing 
on asymptotic SNR regimes.) Throughout this subsection, we will denote the spatial correlation matrix, 
either receive or transmit, by the generic symbol L. 

Before addressing the achievable sum rate of MMSE receivers, it is convenient to give the following 
new result for the ergodic mutual information with optimal MIMO receivers, which simplifies and unifies 
prior expressions in the literature for semi-correlated Rayleigh fading. 

Lemma 1: Consider the transmit-correlated Rayleigh channel H ~ CA/'Ar^,Art(0, Iat^ ® L), or receive- 
correlated Rayleigh channel H ~ CJ\fNr,Nti^^'^ ® Ijv*)' where the spatial correlation matrix L has 
dimension q x q (i.e. for transmit-correlation q = Nt, for receive-correlation q = N^), with eigenvalues 
Pi > ■ ■ ■ > Pq. Also, let p E {Nr, Nt}\q. Then the ergodic MIMO mutual information of H with isotropic 

^Note that in some cases deviations from tliis model have also been observed [41]. 
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inputs and optimal receivers is given by 



r^\Nr,Nt,sm)= y J2 detEp,,(A;) (59) 

where Ep g(k) is a g x g matrix with (s,t)th entry 



Proof: See Appendix Hill ■ 
It is important to note that Lemma [7] allows the correlation to occur between the transmit or receive 
antennas, and places no restrictions on the system dimensions. This is in contrast to prior analyzes (see 
[8, 9, 45]) which have given a separate treatment depending on whether the correlation occurs at the end 
of the link with the least or most number of antennas. 

Given Lemma [7J we can now obtain exact closed- form solutions for the achievable sum rate of MIMO 
MMSE receivers in semi-correlated Rayleigh fading. It is convenient to treat the cases of transmit and 
receive correlation separately. 

Proposition 4: Let H ~ CM j^^^^t^-,'^ ® lAft)^ with L defined as above. Then the MMSE achievable 
sum rate is given by 

LLe<kyyk Pe) \k=Nr-n+l k=Nr-n'+l / 

Proof: The first term in (fTTI) is evaluated directly from Lemma [7] The remaining terms are directly 
inferred from Lemma [7J upon noting that Hj ~ CJ\fNr,Nt-iiOj L ® IjVj_i). ■ 
Proposition 5: Let H ~ CM^^^Nti^^'^Nr ®L), with L defined as above. Also, let L'* denote the (i,i)th 
minor of L, with eigenvalues > ■ ■ ■ > /3i^Nt-i- Then the MMSE achievable sum rate is given by 

N 1 ^* 
I—{Nr,N„snT)= \l^'\ J2 detE^^Mk) 

[h<kyPk - Pe) k=Nt-n+l 

Nt ^ Nt-1 

Nr,Nt-i{k-,i) , (62) 

1=1 [[e<k \Pi,k ~ Pi/) fc=7Vi-?i' 

where ENrM-^ik^^) defined as in (l60l) . but with Pi^k replacing Pk- 

Proof: The first term in (fTTI) is evaluated directly from Lemma [7] The remaining terms are directly 
inferred from Lemma [TJ upon noting that ~ CJ\fNr,Nt-iiO^ '^Nr ® L")- * 
It is important to note that the results in Propositions |?] and \5\ apply for arbitrary numbers of transmit and 
receive antennas. 
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Fig. 3. Achievable sum rate of MIMO MMSE receivers in transmit-correlated Rayleigh fading; comparison of analysis and simulations. 
Results are shown for A'^t = 3 and A^^ ~ 5, and for different correlation coefficients p. 

The result in Proposition \5\ is demonstrated in Fig. [3l where it is compared with the exact MMSE 
achievable sum rate based on Monte-Carlo simulations, for two different transmit-correlation scenarios. 
Here, the simple exponential correlation model was employed, in which case the correlation matrix L 
was constructed with j)th entry p'*^-'', with p denoting the correlation coefficient. We see a precise 
agreement with the analysis and simulations, as expected. Moreover, the MMSE achievable sum rate is 
seen to degrade as the level of transmit correlation is increased, especially in the high SNR regime. 

2) High SNR Analysis: Here we consider channels of the general form (|58l) . The key focus, once again, 
is on the high SNR power offset for the case Nr > Nf. This is given by the following key result: 

Proposition 6: For transmit and receive correlated Rayleigh faded channels, the high SNR power offset 
for a MIMO system with MMSE receiver is given by 

£— (R, S) = £— (liv,,, livj + /(S) + g{R) (63) 
where £J^™^°(Iiv^, ItvJ is the power offset in the absence of spatial correlation given in (|4TI) . and /(■) and 
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Fig. 4. Shift in the high SNR power offset of MIMO MMSE receivers due to spatial coiTelation. Results are shown for both transmit and 
receive correlation, as a function of correlation coefficient p, for different Nr — Nt = n. 



g{-) are given by 



and 



1 ^* 
* k=l 



k,k 



^^ 1 det Yjv,._jv,+i(r) 



(64) 



(65) 



respectively. Here, r = (ri, . . . ^r^r)'^, with ri > . . . > rj^^, are the eigenvalues of R, and YAr^_7Vt+i(r 
denotes an A^^ x A^^ matrix with {s,t)th element 

r*-i for tj^ Nr- Nt + 1 

r*-^ log2 rs hi t = Nr- Nt + 1 



Nr-Nt+l{'''^))s,t 



(66) 



Proof: The result is easily obtained from (I201) . upon invoking the following result [46] 

Nt 

En [J{Nr, Nt, H)] = log2 det S + log2 ej^i^i^) 



and noting that ~ CU Nr,Nt-^{'^^ ^ ® S^^), where S'''' is the (fc, A;)th minor of S. 



(67) 



Thus, in the high SNR regime, the effect of both transmit and receive correlation are clearly decoupled, 
being determined by the functions /(•) and g{-) respectively. Considering the case of transmit correlation, 

*Note that an equivalent expression for ([67} can be found in [38]. 
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we see that 

/(S)>-^log2det-i(S)>0 (68) 

with equality for S = Ia^^. This result can be obtained from (l64l) upon noting that (i) the set of diagonal 
elements [8^^]^ are majorized by the set of eigenvalues of S^^, which, in turn are the reciprocals of the 
eigenvalues of S, and (ii) the function ^ . log2 Xi is Schur-concave. Thus, we see from (l68l) that transmit 
correlation, whilst not affecting the high SNR slope (fT9]) . reduces the achievable sum rate of MMSE 
receivers in the high SNR regime through an increased high SNR power offset; as already observed 
experimentally in Fig. |3l The high SNR approximation based on (fT3l) and (|63l) is also presented in Fig. 
|3l and is seen to converge to the exact MMSE achievable sum rate for quite moderate SNR levels. 

In contrast to the case of transmit correlation, the overall impact of receive correlation is not immediately 
evident from g(-) in (|65l) . due mainly to the presence of the Vandermonde determinant in the denominator 
of the second term. This expression does reveal, however, that the effect of R is purely through its 
eigenvalues, and, interestingly, the relative impact of receive correlation not only depends on Nr, but also 
on Nt. This is in contrast to the effect of transmit correlation in /(•), which depends only on Nt. 

Fig. m plots the shift in high SNR power offset (in dB) due to transmit correlation, based on (|64l) . and 
the shift due to receive correlation, based on (l65l) . as a function of the correlation coefficient p, where S 
and R are constructed according to the exponential correlation model with (i, j)th elements Sij = p'*"-'! 
and Rj j = p'*"-'' respectively. From the figure, we can conclude that the SNR penalty for a MIMO MMSE 
system increases with the level of transmit or receive correlation. Interestingly, the figure also shows that 
for a given correlation coefficient p, the SNR penalty (for the Nr = Nt = n scenarios considered) is more 
severe if the correlation occurs at the transmitter, rather than the receiver; with this difference being most 
significant for small n. 

Now consider the excess high SNR power offset. To evaluate this, we require the high SNR power 
offset with optimal receivers which for the transmit-receive correlated case was first presented in 
[38, Eq. 28]. Using a result from [46], an alternative simplified expression can be obtained, as given by 
the following lemma. 

Lemma 2: For transmit and receive correlated Rayleigh faded MIMO channels, the high SNR power 
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offset with optimal receivers is given by 

which, for the special case A^,. = Nt = n, reduces tC 

= log2 n - log2 (log2 det S + log^ det R) . (70) 



The excess high SNR power offset is now readily obtained from (1691) and (|631) . 

Corollary 8: For transmit and receive correlated Rayleigh faded MIMO channels, the excess high SNR 
power offset is given by 

A^, - 1 

Aex = - log2 e^-— + (71 (S) + (72 (R) (71) 

where 

^1 (S) = ]^ (^E [S-^] + log, det S j (72) 

and 

Ef:«,-A.,+2 det Y,(r) - (iV, - 1) det YA-„_«.+,(r) 
= N.n,.Jr,-r,) • 

From (|68l) and (|64|) . it is easy to establish that (|72|) is non-negative, ie. (71(8) > 0, indicating that in the 
high SNR regime MMSE receivers incur a more significant rate loss due to transmit correlation, compared 
with optimal MIMO receivers. 

3) Low SNR Analysis: For our low SNR analysis, we consider channels of the general form (|58l) . In 
this regime, the main focus, once again, is on characterizing the wideband slope. 

Proposition 7: For transmit and receive correlated Rayleigh faded channels, the wideband slope for 
MIMO MMSE receivers is given by 

^mmse 2,Nj-Ni (74) 

" (2iV, - l)C(R) + (iV,C(S) - ^ E£i C(S"0) 
where S** is the (i,i)t\\ minor of S. 

Proof: For the correlated Rayleigh fading model (|58l) . we can infer the following from [7], 

C(HHt) = C(R) + ^C(S) , 

'This special case was also reported in [38, Eq. 28]. 
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C(H,Hl) = C(R) + ]v^C(S^O • (75) 

Substituting dVS]) into (l33]) leads to the result. ■ 
Note that since the diagonal elements of both R and S are unity, it follows from (|23l) that the dispersion 
numbers in (1751) particularize to 

C(R)-^^, C(S)-^, C(S)-^^^. (76) 
For the case of receive correlation only (i.e. S = IatJ, (l74l) admits the very simple form 

2NrNt 

{2Nt - l)C(R) + N, 

Since C(R) satisfies 



-^o"""'" = , . (77) 



1 < C(R) < A^r (78) 

with the lower bound achieved if the antennas are uncorrelated and the upper bound achieved if the 
antennas are fully correlated, we see from ([771) that receive correlation reduces the achievable sum rate 
of MMSE receivers in the low SNR regime, as quantified by a reduction in wideband slope. It is also 
interesting to compare (1771) with the wideband slope for MIMO with optimal receivers, given by [7] 



Thus, we have the ratio 

5"— _ NtC{R) + Nr 



S'°P* {2Nt - l)C(R) + Nr 
which, based on (fTSl) . also decreases with receive correlation, satisfying 

Nt + l ^ S^^'^ ^ Nt + N, 



(80) 



< ^ < w^r^^r^- (81) 





This result indicates that not only are MMSE receivers degraded at low SNR due to receive correlation, 
but they actually incur more of a loss than do optimal MIMO receivers. Interestingly, we also see that 

2 

and, for C(R) 0, 



^mmsc 



lini ^ = 1 (82) 



Cmmse i 

lim ^ = ^ , (83) 



which is the same limiting behavior observed previously for uncorrelated Rayleigh channels in q5 
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For the case of transmit correlation only (i.e. R = Iat^), focusing on the scenario Nt = 2, (|74l ) reduces 

to 

47V 

° " 3 + iV,(2C(S)-l) ■ ^ 
This result reveals that in the low SNR regime, the effect of transmit correlation in the channel mirrors 
that of receive correlation by reducing the MMSE achievable sum rate through a reduction in wideband 
slope. 

C. Uncorrelated Rician Fading 

We now particularize the general results of Section |lll] to Rician fading channels, representative of line- 
of-sight environments. For convenience, we focus on uncorrelated Rician channels with rank-1 specular 
component, described by 



H ~ CUn.^n. I \l -^a(e,,)a^(^*), J^^n. ® 1 (85) 

where K is the Rician i^'-factor, and a(-) denotes an array response (column) vector (see [47, Eq. 5]), 
parameterized by the angle of arrival Oj. and angle of departure Ot of the specular component. 

It is important to note that, to our knowledge, the results in this section present the first analytical 
investigation of the achievable sum rate of MIMO MMSE receivers in the presence of Rician fading. 

1) Exact Analysis: An exact expression for the MMSE achievable sum rate can be easily obtained by 
evaluating (fTTI) using exact results for the ergodic mutual information of MIMO Rician channels with 
optimal receivers, given in [10,39,44]. We choose to omit explicitly presenting this result here to avoid 
the introduction of more cumbersome notation. 

2) High SNR Analysis: In Rician fading, the high SNR power offset (for the case A^^ > Nt) is given 
by the following key result: 

Proposition 8: For uncorrelated Rician faded channels, the high SNR power offset for a MIMO system 
with MMSE receiver is given by 

£™(ir) = /:~(0) + h{K) (86) 
where /:™(0) is given by (gB and 

h{K) = \og,{K+l) - K\og^e{Nte{Nr,Nt,K) - {Nt - l)e{Nr,Nt - 1,K)) , (87) 
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15 

Rician K-Factor 



Fig. 5. Shift in liigh SNR power offset hi{K) and excess power offset h2{K) in Rician fading. Results are siiown for A''^ = Nt = 2. 



with 



e{Nr, Nt, K) = 2^2(1, 1; 2, Nr + 1; -KN^Nt, 



(88) 



where 2-^2 (■) denotes the generalized hypergeometric function. 

Proof: The result is easily obtained from (|20l ). upon invoking the following result [38] 

Nt-l 

En [JiN,, Nt, H] = log2 e ^(^^ - i) - N^ \og,{K + 1) + KN^ log, ee{Nr, iV„ K) (89) 
and noting that 

Hfe ~ CN 



1=0 



K 



1 



-iNr ® IjVt-1 I ; 



(90) 



K+l ' ' " K + 1 

where afc(-) corresponds to the response vector a(-) with the kth element removed. ■ 
Interestingly, we see that the impact of line-of-sight in the high SNR regime is only through the relative 
strength of the specular component (ie. through the Rician i^-factor), and is independent of its geometry 
(ie. independent of 9r and 9t). Moreover, Proposition\^re\eals that line-of-sight imposes a shift in power 
offset, as succinctly characterized through the function hi{K). This behavior is investigated in Fig. [51 
where we plot hi{K) (in dB) as a function of K. We clearly see that the high SNR power offset increases 
monotonically with K, revealing that the presence of line-of-sight reduces the achievable sum rate of 



26 

MIMO MMSE receivers in the high SNR regime. 

Corollary 9: For uncorrelated Rician faded MIMO channels, the excess high SNR power offset is given 

by 

AexW = Aex(0) + /l2W (91) 

with Aex(O) given by (gV]), and 

h2{K) = - log2 eK{Nt - 1) {e{Nr, Nt, K) - e{Nr, Nt-1,K)) . (92) 



Proof: This result is obtained by substituting (1861) and [38, Eq. (67)] into (|22l) . ■ 
As shown in Fig. [51 the function h2{K) increases with K, confirming that not only does line-of-sight 

decrease the high SNR achievable sum rate of MMSE receivers through an increased power offset, but 

the loss is more significant compared with optimal MIMO receivers. 

3) Low SNR Analysis: For Rician channels in the low SNR regime, we have the following key result: 
Proposition 9: For uncorrelated Rician faded channels, the wideband slope for MIMO MMSE receivers 

is given by 

° K^{2Nt-l)Nr+i2K+l){2Nt + Nr-l)' 

Proof: For the uncorrelated Rician fading model (|85l) . we can infer the following from [7], 

j^^j^2 _^ {Nr+Nt){2K+l) 



C(HHt) 



Nt 



(K + l) 



-2 I {Nr+Nt-l){2K+l) 



C (h-hI) = ■ <''^) 

Substituting ^ into (|33]) leads to the result. ■ 
As also observed previously for the high SNR regime, we see that the impact of line-of-sight in the 
low SNR regime is only through the relative strength of the specular component (ie. through the Rician 
iT-factor). Moreover, the wideband slope (l93l) is a monotonically decreasing function of K, implying that 
line-of-sight has a damaging effect on the achievable sum rate of MIMO MMSE receivers in the low SNR 
regime. Comparing this result with the corresponding wideband slope for MIMO with optimal receivers, 
given for Rician fading by [7] 



27 



we obtain the interesting relationship 



(96) 



5°^* cp{K,Nr,2Nt-l) 



with 



ip{K, m, n) = K'^mn + {2K + l)(m + n) . 



(97) 



This ratio, again, is a decreasing function of K, satisfying 



(98) 



where the lower bound is approached as X — > oo, corresponding to the purely deterministic channel 
scenario, and the upper bound is approached as — > 0, corresponding to Rayleigh fading. 



We presented a new analytic framework for investigating the achievable sum rate of MEMO systems 
employing MMSE receivers, revealing a simple but powerful connection with the ergodic MIMO mutual 
information achieved with optimal receivers. This framework allowed us to directly exploit existing MIMO 
results in the literature, thereby circumventing the major challenges entailed with explicitly characterizing 
the SINR distribution at the MMSE output. To demonstrate the utility of the framework, we presented 
particularizations for uncorrelated and correlated Rayleigh fading, and uncorrelated Rician fading channels, 
yielding new exact closed-form expressions for the MMSE achievable sum rate as well as simplified 
expressions for the high and low SNR regimes. Through these expressions, we obtained key analytical 
insights into the effect of the various system and channel parameters under practical fading conditions. 
For example, we demonstrated that at both high and low SNR, the MMSE achievable sum rate is reduced 
by either spatial correlation or line-of-sight. At high SNR, this rate reduction is manifested as an increased 
power offset, whereas at low SNR, through a reduced wideband slope. Moreover, at both high and low 
SNRs, the rate loss due to spatial correlation or line-of-sight was shown to be more significant for MMSE 
receivers than for optimal receivers. We also demonstrated that the effect of line-of-sight on the MMSE 
achievable sum rate was dependent on the relative strength of the specular component, but not the geometry 
of such component. 

We would like to stress that although the main focus of the paper was on single-user MIMO systems 
with MMSE receivers, many of the results apply almost verbatim to multi-user scenarios; in particular. 



V. Concluding Remarks 
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the analysis of multiple access channels with MMSE receivers, and MIMO broadcast channels with 
either MMSE-based transmit precoding (see, eg. [23]) or MMSE reception [48]. Moreover, the proposed 
framework extends to many other scenarios beyond those explicitly studied in this paper. These include, 
for example, single-user MIMO systems operating in the presence of interference [19, 20, 22] and amplify- 
and-forward relaying systems [46]. 

Appendix I 
Proof of Theorem [I] 

We start by substituting ^ into ^TQ, and usin^ [50] 

detZ 



' det Z* 

where Z** is the (z,i)th minor of the matrix Z, to yield 



(99) 



Nt 



snr 



i=l 



log2 



det ( Ijv, + f;^HtH 



det + 



t-C/H 



logo det {lm + ^H"^H 



i=l 



log2det(l^,_i + ^(HtH)' 



(100) 



Noting that (H+H)** = H^Hi, the resuk follows from ©. 



Appendix II 
Proof of Theorem [3] 



From (fTTI) . and noting that 



we can obtain 



— In det (I + xA) 

dx 

In det (I + X Al 



d^x 



x=0 



x=0 



tr (A) , 
= -tr (A^) 



Nt 



/~(0) = log2 e ( Eh [tr (HRt)] - - E^^ 



tr H,H 



(101) 



(102) 



Note that this property has also been used in relation to linear MIMO receivers in [49]. 



29 



and 



Nt 



iV,EH[tr([HHt]2)]-5^EHJtr [H,Hl 



Due to the channel power nomialization we have 



i=l 



En [tr (HHt)] = NrN^, En, 



2 Nt 



i=l 



tr ( HiHl 



iV,(iVi - 1) 



regardless of the specific channel statistics. As such, (11021) evaluates to 



In 2 



(103) 



(104) 



(105) 



which, from ([28]), yields (|32l). The result ^ is obtained by substituting (11051) and (11031) into ^ and 
simplifying. 



Appendix III 
Proof of Lemma [H 

The MIMO mutual information with isotropic inputs is given by 

P''\Nr,Nt,sm) = En 
= n 



logadet f + — HtH 
/°^log2(l + ^A)/.(A)dA 



(106) 



where A is an unordered non-zero eigenvalue of H^^H, with probability density function (p.d.f.) /a(-)- 
Recently, the following unified expression (applying for arbitrary and Nt) was presented for this p.d.f. 
[46] 

1 



/a(A) 

where is a g x g matrix with entries 

{D = 



det Da 



(107) 



t ^ k, 

x^Zl^p-^/Psm-P-i f - I. 
The result follows by substituting (11071) into (11061) and integrating using an identity from [10]. 



(108) 
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